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Abstract: We present an efficient, covariant, graph-based method to integrate superfields 
over fermionic spaces of high dimensionality. We illustrate this method with the compu- 
tation of the most general sixteen-dimensional Majorana-Weyl integral in ten dimensions. 
Our method has applications to the construction of higher-derivative supergravity actions 
as well as the computation of string and membrane vertex operator correlators. 
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1 Introduction 

Despite the conceptual elegance of superspace methods, their use in the study of supergravity and 
string theory is hampered by a number of technical difficulties. The problem which is perhaps most 
manifest is the fact that high-dimensional fermionic integrals, although conceptually simple, are hard 
to evaluate explicitly. This forms a serious obstacle when one attempts to relate the often elegant 
superspace expressions to results in terms of supergravity component fields. 

The fermionic integration problem is encountered in a variety of different situations. One of these 
is the construction of ten-dimensional higher-derivative supergravity actions. By virtue of the pres- 
ence of the dilaton field, these supergravity theories allow for the existence of scalar superfields which 
contain the entire supergravity multiplet in their component expansion 0,1^. This observation has 
led (now almost twenty years ago) to the hope that complicated higher-derivative actions can perhaps 
be constructed in terms of superspace integrals of simple expressions |^, Q|. Although the application 
of this idea to the type-IIB theory is beset with difficulties [5|, even the simpler construction in = 1 
supergravity has never been worked out in full detail. One particular problem which has remained 
unsolved is how to relate superspace expressions to component ones. This same problem also appears 
in the computation of vertex operator correlators for superparticles 0|, strings and membranes 0, 
when these are formulated using target-space spinors (i.e. using the Green-Schwarz or Berkovits for- 
malisms). High-dimensional fermionic integrals appear here when one tries to integrate out fermionic 
zero modes of the fields living on the world-volume. 

While generic covariant fermionic integrals are thus so far not known, special cases which exhibit 
additional symmetries are sometimes tractable. One such case is the integral that leads to the sixteen- 
dilatino interaction in the type-IIB theory 0, 0| . The superspace supergravity computation and the 
computation of a sixteen-fermion vertex operator correlator both lead to the trivial integral 

j d'^d {OaX")^'' = eai-aieA'^^ ' • • A-^^^ = 16! • • • . (1.1) 

A slightly more complicated expression is obtained for the four-graviton amplitude in the Green- 
Schwarz formalism, when computed in the light-cone gauge. The resulting fermionic integral can be 
decomposed as the product of two known S0(8) fermionic integrals with the result 0] 

I dW^[{h"'^^){h''^)Rmnrsy = (^8*8 " ^ses) . (1.2) 
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Again, this integral also occurs in an analysis of four-graviton couplings in supergravity. A similar 
result can also be obtained for the four-point correlator of supermembrane vertex operators in the 
light-cone gauge 0]- However, it is clear that (|1.H) and (|1.2() form only the tip of the iceberg. Many 
interesting results wait to be derived once a fully covariant way is established to perform an arbitrary 
high-dimensional fermionic superintegral and express it in terms of Lorentz singlets (Kronecker deltas 
and epsilon tensors). 

Therefore, it is the purpose of this small note to discuss a generic, covariant method for the inte- 
gration of arbitrary functions over high-dimensional fermionic spaces. We demonstrate the feasibility 
of our method by deriving an explicit expression for the most general sixteen-component SO (10) 
fermionic integral (the 'W = 1 integral") in terms of Lorentz singlets. This result is rather interest- 
ing by itself, and we will discuss our motivation to derive it, including possible applications, in some 
more detail towards the end (in section The method can be applied easily to the other ten- and 
eleven-dimensional supergravity theories. Expressions for S0(9) integrals, relevant for superparti- 
cle and supermembrane calculations in eleven-dimensional supergravity in the light-cone gauge, will 
appear shortly [l3 |. 

For completeness, we also describe, in the appendix, an efficient method to reduce tensor polyno- 
mials to a minimal basis. This simple method does not seem to be widely known, but is of considerable 
help in dealing with higher-derivative Lorentz invariants. 



2 Fermionic integrals 

2.1 A simple eight-fermion example 

The goal of this note is to show how high-dimensional fermionic integrals can be evaluated covariantly 
and in full generality. However, the techniques which we will use apply also to much simpler cases. 
It is therefore illustrative to first consider a simpler fermionic integral, which can be done by hand 
and for which the answer has been known for a long time, so as to get familiar with the techniques. 
Let us thus consider the following integral over the eight-dimensional space of S0(8) spinors [lo| . 

jiih-idi j ^^Q± (^iyiii^i^ . . . (g±y4i4^±) ^ (2.1) 

Here the it symbols denote the chirality of the spinors. In order to determine the number of Lorentz 
singlets which is needed to express this integral, we compute the tensor product of the four sym- 
metrised sets of two anti-symmetric vector indices 0] , 

fs X [00001 ©... inS0(8), 
Sym4([010...]) = (B)^,^= (2.2) 

[2 X [00000 . . .] e . . . in S0(2A;) for k > 4. 

This result implies that ()2.1() can be decomposed in two delta singlets and one epsilon singlet (the 
epsilon singlet is dimension dependent and corresponds to the disappearing singlet when the tensor 
product is evaluated in higher dimensions). It is straightforward to find these three independent 
singlets; we will use 

Di = -L (5^1*2 5Jii2<5»3i4ji3i4 + II ^grms) , 

D2 = ^ (5iii2 ji2i3^i3i4 ji4ii + 47 terms) , (^-^^ 
= ^nji«2---i4 _ 

The fermionic integral can thus be written as 
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and the goal is to determine the unknown coefficients ai,a2 and /3. 

By using an exphcit representation for the S0(8) gamma matrices, it is straightforward to evalu- 
ate H'2.1() for particular values of the eight indices (by selecting the terms in the resulting polynomial 
of spinor components in which each component occurs once). Similarly, it is straightforward to deter- 
mine the value of Di, D2 and E for a particular set of index values. Three independent combinations 
are hsted below, 



[iiji]--- [Uj4] 


/+ 


/_ 


Di 


D2 
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[12] [12] [34] [34] 


-128 


-128 


1/12 








[12] [23] [34] [41] 


128 


128 





1/48 





[12] [34] [56] [78] 


128 


-128 








1 



(2.5) 



This leads to three equations for three unknowns, from which one determines the coefficients to be 

ai = -1536, 02 = 6144, (3 = 128. (2.6) 
Comparing with the tg tensor of appendix 9. A of ^3]) one then obtains the expected result 

/± = 256i8,±. (2.7) 

In the next section we will see that the sixteen fermion integral can be evaluated using precisely 
the same logic, although the number of Lorentz singlets increases sharply and it also becomes more 
complicated to evaluate their values given a set of indices. This increased complexity calls for a 
number of new ideas. 



2.2 The sixteen fermion integral 

Let us now turn to the evaluation of the sixteen fermion integral in ten-dimensional simple super- 
gravity. The Weyl and Majorana properties of the spinor 9 imply that bilinears in 9 can be written 
as three- forms. The most general integrand therefore has the form 

..|,..,(5r...^.,)(5r.«*,)...(,-r.«.'..). ,2.8) 

The goal will again be to express this integral in terms of Lorentz singlets, i.e. Kronecker deltas and 
epsilon tensors carrying the free vector indices. That is, we want to write the integral as 

i 

in terms of a set of basis tensors T(j) . 

The first step in this integration is to determine the number of Lorentz singlets in which the 
integral ()2.8|1 can be decomposed. This number is easily obtained by considering the tensor product 

/ x8 f 33 X [00000] ©... inSO(lO), 
Sym8([00100...]) = (|) = (2-10) 
V /sym 1^24 X [000 . . .] ® . . . in S0(2A;) for A; > 5. 

There are thus 33 singlets in SO(IO). Nine of these disappear when one considers higher-dimensional 
spaces, and these thus correspond to parity-odd invariants involving the ten-dimensional epsilon 
tensor. 

In contrast to the situation in the previous section, it is now not so easy to guess the 33 Lorentz 
singlets. However, the explicit construction of these singlets can be translated to an elegant problem in 
graph theory. Let us first focus on the singlets which involve only Kronecker deltas. We will represent 
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fully anti-symmetrised index triplets by trivalent nodes, and represent Kronecker deltas which set two 
indices equal by edges connecting the nodes. Multiple edges are allowed (e.g. S^^^2§jij2§kik2 ^ which 
corresponds to two nodes with a 3- fold connection between them). The problem of finding all 24 
independent singlets now corresponds, in standard graph terminology, to the problem of finding all 
3-regular multigraphs (not necessarily connected) with eight vertices. All graphs of this type can be 
found with the help of 3; ^i and one finds a total of 32 graphs. This may seem to contradict (|2.10)) . 
however, the anti-symmetry of the index triplets is not yet fully encoded at this stage. Using anti- 
symmetry, 8 of the 32 graphs can be shown to correspond to a vanishing expression. For example, 
the double contraction of two identical three-forms is symmetric in the two free indices and therefore 
vanishes when contracted into a further three-form. In graphical notation, this is seen from 







0. 



(2.11) 



Another vanishing subgraph is 






(2.12) 

These identities make 8 graphs vanish identically, see figure |2I The remaining 24 graphs are listed in 
figure Hand their explicit expressions in terms of Kronecker deltas can be found in tabled 

The construction of the graphs corresponding to parity-odd singlets is less systematic. We can 
again introduce a graphical notation, representing an epsilon tensor by a 10-valent vertex (or "box"). 
Many graphs can be formed out of one 10-valent and eight trivalent vertices, but nine of these are 
sufficient to represent the nine independent parity-odd singlets in the tensor product (12.10(1 . These 
are displayed in figure El and the corresponding explicit expressions can be found in table |2j 

Having determined the basis on which (|2.8|) can be decomposed, the remaining step is to determine 
the coefficients in front of each of these basis tensors. This is again done by matching the values of (|2.8() 
and Ij2.9j) for various sets of values of the indices, and solving the resulting system of linear equations.^ 
In graph-theory language, the evaluation of (|2.9|) corresponds to finding all ways of colouring the edges 
of a graph with numbers, given a set of three numbers at each vertex. This requires some care in 
order to keep the computation within bounds. The basis tensors in table ^ and [21 contain an implicit 
anti-symmetrisation over all indices in each triplet, as well as an implicit symmetrisation over all 
index triplets. A brute- force algorithm which investigates the value of a given singlet for all terms in 
the symmetrisation therefore leads to a worst-case situation in which (3!)^ • 8! ~ 6.8 x lO^'' different 
terms (or colourings) have to be considered. The graphical representation suggests a much more 
efficient "backtracking" algorithm. This algorithm constructs the graph labellings vertex by vertex 
and checks after each choice if the index assignment is still consistent. If it is, the next vertex is 
labelled, if not, the algorithm backtracks and proceeds to the next choice of edge labelling. 

The result of this matching and the subsequent solution of the linear system is presented in the 
last column of tables ^ and [3 This concludes the computation of ()2.8|) . 



2.3 Covariant computation of the i?^ integral 

A useful check of our method is to compute the well-known term in the linearised heterotic or 
type-IIB theory, and verify that it reproduces results which were obtained previously using non- 



^The chirality of the spinor 8 in 12.81 enters at this stage. As in section ITTI the choice of chirahty reflects itself in 
the sign of the coefficients 025 . . . 033 of the parity-odd singlets, and we will not comment on this any further. 



5 



Graph i 
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Table 1: The cocfEcients of the 24 Lorentz singlets in the 0^^ integral which can be expressed solely 
using Kronecker deltas. The singlets are understood to be antisymmetrised in the [ijk] triplets, and 
symmetrised over 1 ... 8. 
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Figure 1: Graphical display of the 24 Lorentz singlets in the 9^^ integral which can be ex- 
pressed solely using Kronecker deltas. Each dot represents a fully anti-symmetrised index triplet 
(e.g. and the lines indicate how these indices appear on Kronecker delta symbols. 



covariant methods. The i?'* term arises from a superfield integral of the form 
where the relevant terms in the scalar superfield $ take the form 



(2.13) 



Hjfct -r • ■ . . (2-14) 

By making explicit use of the results of the previous subsection, we can express the integral in terms 
of the 24 independent parity-even Lorentz singlets, 

24 

niiiifci-igjsfcg [„. . p. . . . 1 . . . . p. , . . 1 (2.15) 

i=l 



(Note that the parity-odd part of the integral does not contribute, since there are no ei?^ scalars 
in ten dimensions.) The individual terms in the sum lead to lengthy linear combinations of the 
26 independent quartic curvature scalars. All dependence on the Ricci curvature cancels in the total 
result, and by decomposing the result on the basis of the 7 Fulling invariants (see also appendix^ 
we are left with the required result in terms of Weyl tensors, 

SlOSlo)^^^ — 2^"^?,^ { — \C^'^^^Cpq^Crt'^Csuvw + C'^'^'^^ Cp^ Ct' Cuvsw) ■ (2.16) 

Imposing this result actually turns out to be restrictive enough to completely determine the coefficients 

ai . . . Q!24- 
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Figure 2: Graphs of parity-even type which vanish identicaUy by virtue of the identities l|2.11|l 
and (EEl). 



3 Discussion and applications 

We have presented the computation of the ten-dimensional fermionic integral (|2.8j) in terms of 33 
basis tensors, using a general method applicable to all fermionic integrals, in particular those of ten- 
and eleven-dimensional supergravity. The solution of this problem was facilitated by mapping it onto 
a graph construction and colouring problem. To conclude, let us discuss a number of applications of 
this integration procedure. 

Let us start with the = 1 heterotic theory. The use of the dilaton superfield for the construction 
of the i?^ invariant was proposed a long time ago 01 • Although the expansion of the scalar superfield 
was only worked out to lowest order, it should be possible to extend this analysis to higher order in 9 
using computer assistance (a similar analysis would have to be performed to construct the measure). 
One particular reason why it is interesting to pursue this program is that it could provide insight into 
the still elusive supersymmetrisation of the {tsts + |eio£io)-R^ invariant. As was emphasised in [l^ . 
the tgts part of this construction is reasonably well understood because of its relation to a super- 
Yang-Mills invariant. The eeR'^ term in contrast, does not admit a derivation by "squaring" a 
super- Yang-Mills action. However, this term potentially plays an important role in the modification 
of the superspace torsion constraints. A direct derivation in components would certainly help to 
understand this issue. 

The complete set of higher derivative interactions that accompany the R'^ term involve, in par- 
ticular, fluxes that are essential for understanding nontrivial compactifications of string/M-theory. 
Unfortunately, there are intrinsic difficulties in the superspace description of such interactions due to 
the absence of an off-shell superspace formalism for theories with maximal supersymmetry. However, 
a limited amount can be deduced from the available on-shell superspace formulations. For example, in 
the type-HB theory the superfield formulation of Q] encapsulates the classical theory which contains 
component fields comprising the one- form (made from a derivative of the complex scalar), the NS-NS 
and R-R three- forms and the five- form field strength (as well as the fermions). There are difficulties 
in extending this to the general leading higher derivative interactions [3| but in the special case in 
which only the five-form and the metric are non-vanishing, the five-form dependence enters purely as 
a companion to the curvature in the 6^ term in the scalar superfield 0, ^y]. As shown in this 
leads to an elegant understanding of the nonrenormalisation of the Z)3-brane supergravity solution 
by these leading higher derivative interactions. In other examples string perturbation theory has 
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Figure 3: Graphical display of the 9 Lorcntz singlets in the 9^^ integral which contain epsilon 
tensors. Rectangles denote epsilon tensors, and each dot represents a fully anti-symmetrised index 
triplet, as in figure ^ 
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Vhii 
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Table 2: Coefficients of the nine selected epsilon singlets which occur in the fcrmionic integral (|2.8|l . 
The vanishing of four of the coefficients has been achieved by choosing a suitable basis of the nine 
parity-odd singlets in (|2.10|l . 
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provided evidence for the structure of such terms, but a superfield analysis has not yet been possible. 
For example, certain H^E? terms were determined from a string calculation in 2^, 3. In any case. 



the technical results of this paper may be of value in any future progress towards an understanding 
of the superspace formulation of such interactions. 

A completely different application of our superintegration techniques concerns the study of vertex 
operator correlators. In space-time super symmetric formalisms for superparticles, strings and mem- 
branes, the leading higher-derivative amplitudes arise as simple expressions over world-sheet fields, 
involving only a fermionic zero-mode integral. Such integrals are precisely of the type considered here. 
Using the methods of the present paper, it has recently become possible to determine the {DFu^E? 
and (L'F(4))^ terms in the M-theory effective action directly from a superparticle calculation 
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A Appendix: reduction of tensor polynomials 



Once one has done a superspace integral over a sufficiently complex integrand, one typically ends 
up with tensor polynomials which can be further reduced using the symmetries of the individual 
tensors, as well as the exchange of identical tensors. We would here like to comment briefly on a 
reduction method which incorporates all symmetries, including the Ricci cyclic identity and Bianchi 

0, 



identities (this method was used implicitly in several papers, see e.g. |22l. 121^ . but as far as we know 
was never spelled out, and more importantly, does not seem to be widely known; it is not mentioned 
in the standard reference and is also missed in most recent literature on tensor polynomial 
simplification) . 

The simplest symmetries of tensors are mono-term symmetries, such as anti-symmetry in a set of 
indices. These symmetries always relate one particular index distribution to one other distribution, 
and can be used step-by-step to reduce a given tensor monomial to a canonical form. Various efficient 
algorithms have been discussed in the literature and we will not comment on these symmetries any 
further (see 2^, 2^, 2^ and in particular the implementation in H^]). 

The more complicated symmetries are multi-term symmetries, such as the Ricci cyclic identity or 
the Bianchi identity, which relate a sum of terms with different index distributions. These symmetries 
are all manifestations of the so-called Garnir symmetries of Young tableaux |2^. These state that 
the sum over all anti-symmetrisation of boxes in a Garnir hook is identically zero. Examples of such 
Garnir hooks are given below. 



which represent the Ricci cyclic identity, the Bianchi identity on a five-form and a more general Garnir 
symmetry, respectively. Applying a Garnir symmetry on a tensor produces a sum of tensors, which 
means that one can no longer restrict to the canonicalisation of tensor monomials. 

There is, however, a simple way to construct instead a basis of monomials which takes the Garnir 
symmetries into account, and reduce any given expression to this basis. It consists of simply replacing 
each tensor in a monomial by its Young-projected form. A basis of monomials can now be constructed 
by first generating a list of all possible index contractions of the tensors, and then projecting each of 
these using the Young projection on each of the individual tensors. Monomials which are identical 
by virtue of Garnir symmetries will then map to the same sum of monomials. 

There are, however, a few subtleties. Firstly, it would be prohibitively expensive to write down all 
terms in a Young-projected tensor. Instead, it is much more efficient to reduce the Young-projected 
forms bv making use of the mono-term symmetries, which are easy to deal with using the methods 
of 0,0,0. One thus obtains, for e.g. the Riemann tensor, 



^abcd 



2 (2 Rabcd 



Radbc + Racbd) , 



(A.l) 



instead of the (2!)^ = 16 terms which are produced by the Young projector. The expression 
on the right-hand side manifestly satisfies the cyclic Ricci identity, even if one only knows about 
the mono-term symmetries of the Riemann tensor. Using the projector (|A.1|) it is easy to show 
e.g. that 2RabcdRacbd = RabcdRabcd- The mouomial on the left-hand side maps to 



RabcdRacbd ^ 3 (^RabcdRacbd ~l~ RabcdRabcd^ > 



(A.2) 



and similarly RabcdRabcd maps to twice this expression, thereby proving the identity in a way which 
easily extends to much more complicated cases. 
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Secondly, the Young projectors are not dimension dependent. As a result, not all relations between 
monomials will be recognised. The additional relations exist because invariants sometimes arise as 
the contraction of tensors with two d-dimensional epsilon tensors. When written out in terms of 
Kronecker deltas, this leads to anti-symmetrisation in d indices when none of the indices of the 
epsilon tensors are contracted with each other. Clearly, such an invariant can be written down in 
lower dimensions, but it vanishes identically, although this is not recognised by simply performing 
the Young projections. Relations obtained in this way have to be taken into account separately, but 
are easy to find. 
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